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Abstract 

We investigate, at both zero and finite temperature, the properties of strangelets 
versus the electric charge Z and strangeness S. The strangelet radius is not a monotonic 
function of either charge or strangeness, and a minimum is reached in the [Z, S) 
plane. However, the thermodynamically stable strangelets do not correspond to the 
radius minimum. The minimum radius always appears at positive strangeness, while 
the stable radius may appear at negative strangeness for very small baryon numbers. 
For large baryon numbers, the stable radius is proportional to the cubic root of baryon 
numbers, but inversely proportional to the square root of the confinement parameter in 
the present model. If bulk strange quark matter is absolutely stable, the reduced size 
of strangelets is about 1 fm, which may be relevant for the analysis of the strangelet 
production and detection. 



* E-mail address: gxpeng@ihep.ac.cn 



1 Introduction 



Strange quark matter (SQM) has been one of the hot topics in nuclear physics since Witten's 
conjecture pQ that SQM might be the true ground state of strong interactions. Lumps of 
SQM are customarily called strangelets [2J, or for short, slets [3]. There are basically two 
kinds of slets: one is the ordinary slets without pairing [21 El EE El E], the other is the 
color superconductivity slets 0, [3], [10] • The possible terrestrial production of slets has been 
studied in high energy heavy ion experiments, e.g., at the CERN SPS energies For 
a review on recent experimental searches of quark-gluon plasma at Brookhaven RHIC, see 
Ref. [12]. At high densities, quark matter may be the most favorite phase, and could thus 
exist in the core of compact stars [13] . Strange star collisions could release slets as part of 
the energetic cosmic rays, and some of the cosmic ray slets might be on the way to our Earth 
[T4"[ [T5] . The Alpha Magnetic Spectrometer (AMS-02) [16], which is planned to operate on 
the International Space Station, may offer the opportunity to detect these slets |17j . 

It is still an interesting open problem whether or not these cosmic ray slets incident on 
top of the Earth's atmosphere can reach to the ground or sea level. In fact, several exotic 
cosmic ray events with anomalously low charge-to-mass ratio have been observed at rather 
low altitudes [T8l [191 [20] . In literature, one finds a number of possible scenarios of slet 
propagation which depends strongly on assumed size of slets. For example, the radius of 
slets used in [21] is different from that in [221 122] by nearly an order of magnitude. Since 
the mean free path of a slet in the atmosphere is strongly dependent on its radius i? s iet, 
as, for instance, in the simple law A = A a i r mN/[vr(i? s i ct + -R a ir) 2 ] where m-^ rs 939 MeV is 
the nucleon mass, A ah « 14.5 and R ah = 1.12 A 1 ^ « 2.73 fm are, respectively, the mean 
mass number and radius of the nucleus in the atmosphere, different radii led to significantly 
different conclusions. 

Meanwhile, the charge property of slets is very important, as shown by Madsen et al. 
who found that the slets of low charge-to-mass ratio are favored in the ultra-high-energy 



2 



cosmic rays [24J, and the color flavor locked (CFL) slets have charge Z w 0.3A 1//3 where A is 
the baryon number. Jaffe et al. also demonstrated that the slets with large baryon number 
have positive charges Z ~ A 1 / 3 [I]. In addition to these positively charged slets, negative 
charges are also possible for both ordinary [25] and CFL slets [3] in beta equilibrium. At 
the same time, the strangeness fraction is also an important factor to determine the stable 
configuration of slets [6]. 

Recently, we suggested a new quark mass scaling [26] based on the linear confinement, 
and the new model was applied to investigating the properties of slets in full beta equilibrium. 
In heavy ion collision experiments, however, the time scale is not enough for perfect beta 
equilibrium. Moreover, the charge/strangeness composition, and especially the size of slets, 
are quite important and useful to analyse possibility of production and detection of slets 
|21j . In this paper, therefore, we study the relevant properties of slets with the new quark 
mass scaling in [26], without imposing beta equilibrium, and with focus on the slet size. 
It is found that the mechanically stable radius of a slet with fixed baryon number and 
temperature is not a monotonic function of either charge Z or strangeness S. The radius has 
a minimum in the (Z, S) plane. However, the radius minimum does not correspond to the 
minimum of free energy. We determine the composition of thermodynamically stable slets 
by minimizing the free energy. For a large baryon number A, the corresponding slet radius is 
R = [(3/4) 1 / 3 /v / xo-D]^4 1 ^ 3 , where xq = 1.3278478, and D is the confinement parameter in the 
present model. If SQM is absolutely stable, the reduced stable radius is r s i et = R/A 1 ^ 3 fa 1 
fm. In the conventional bag model, slets always contain strange quarks. In the present 
model, we find that very small slets tend to contain anti-strangeness, and the ratio of charge 
to baryon number increases. 

This paper is organized as follow. In Sec. [2] we introduce the thermodynamic treatment 
with density and temperature dependent quark masses. The properties of strangelets related 
to the strangeness and electric charge at both zero and finite temperature are presented in 
Sec. El A summary is given in the final section HI 
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2 Thermodynamics with density and temperature de- 
pendent quark masses 

We start from the total free-particle thermodynamic potential density 

Q = ^2Q i (T,fi i ,m i ,R), (1) 

i 

where the summation index % goes over u, d, s quark flavors , T is the temperature, m ; and /ij (i 
u, d, s) are the corresponding quark masses and chemical potentials, and R is the slet radius. 
At finite temperature, we treat the anti-quarks as a whole with quarks. The contribution of 
the thermodynamic potential density from the density of state n^(p, mj, R) is given in the 
multi-expansion approach [27] as 



a = —T 



f°° | In [l + e -(vV+"\ 2 ^)/ T 



+ In [l + e -(vV+ m ?+^V T ] I n'^p, m h R)dp, (2) 

where the density of state is 

n>(p, m„ R) = l~ 2 \v 2 - ||arctan + 



3p mi 
-arctan — 



(3) 



2rrii \ p 

The three terms on the right are, respectively, the volume term, surface term [21 H] and 
curvature term [28J. 

To include the confinement interaction between quarks, we treat the quark mass as den- 
sity and temperature dependent, i.e. rrii = m^nx,, T), where = J2i with nj (i = u, d, s) 
being the quark number densities. This means that the mass of quarks and antiquarks varies 
with state parameters in a medium. We can divide the quark mass into two parts: one is the 
current mass m i0 , the other is the interacting term mj, i.e., = m, + mj. In the present 
calculations, we take the quark current masses m U Q = 5 MeV, m^o = 10 MeV and m s o = 
120 MeV, respectively. Because the strong interaction between quarks is a color interaction, 
mi is common for all quark flavors. The key point is how to determine the interaction mj. 



In Ref. [29], the ansatz ra Y = J^[l - (T/T c ) 2 ] was introduced. Because it caused an unrea- 
sonable temperature dependence of the slet radius, another term linear in temperature was 
added [S]. Based on the in-medium chiral condensates and linear confinement, we recently 
derived a new quark mass scaling, which can be expressed as 



D 

mo + — 



nr 



8T 



exp 



T 
-A^ 
T 



(4) 



where A = LambertW(8) ~ 1.60581199632 is a constant, T c = 170 MeV is the critical tem- 
perature. The exponent z was previously taken to be 1 [30]. [3T|. [32] . In order to be consistent 
with the linear confinement, derivations based on the in- medium chiral condensates [33] 
showed that it is more reasonable to take z = 1/3 [3H [35j [36]. The confinement parameter 
D can be constrained to a very narrow range by stability arguments and we take D 1 ^ 2 = 156 
MeV [37]. 

The particle number density for each quark flavor can be derived by the following ex- 
pression 

dtt 



7k 



dfii 



(5) 



The pressure is 



P 



^ R dfl x , dQ drrii 
— — + n b}^ 



(6) 



3 OR i dmi drib 

where the last term is due to the density dependence of quark masses [371 E2]- The partial 
derivatives dmi/dn^ in Eq. (jSJ) can be easily obtained from the quark mass scaling in Eq. (j3J), 
i.e., 
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Accordingly the free energy density of the slets is 

tp n V 9(1 
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At zero temperature, the relevant integrations can be carried out. The quark number 
densities in Eq. <^ become 



vf 9m, 



7T 2 An 2 R 
9m,' 



An 2 R 2 



' K x i + l)arctan(xj) — — X{ + 1 



3tt 



7T 



+ l)arctan(xj) — Xi [ ~^ x i ~ 1 



(9) 



where z/« = y/i? — m? is the Fermi momentum of the quark flavor i and x« = Ui/rrii. The 
free energy density in Eq. (JSl) becomes the energy density 



E 
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And the pressure in Eq. ([6]) becomes 



7i — Xi\Jx 2 + 1 — arcsh(xj) — 2(x 2 + l) 3 / 2 arccot(x ( 
7T + XiJxf + 1 + arcsh(xj) — 2{xj + l) 3 / 2 arccot(xj) 



(10) 



Tfl I 

P = E j^(2x 2 - 3)v/x 2 + l + 3arcsh(xi) - 12z 
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mi 
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7r{3Jxf + 1 - 2) - 2arcsh(xi) - 2(x- + l) 3/2 arccot(x;) 
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(11) 



In the above Eqs. (J9J), (flOl) . and (TTTT) . arcsh(xj) = \n(x + \/x 2 + 1) is the inverse hyperbolic 
sine, arctan(xj) and arccot(xj) are the inverse tangent and cotangent functions. Please note, 
the interaction part of the quark mass scaling, i.e., the second term on the right hand side 
of Eq. (j4j), simply gives mi = D/n^ at zero temperature. 
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3 Properties of strangelets 

Thanks to the pioneer works of Witten and Jaffe et al. [H [2] , we have known a lot about 
slets. The properties of slets away from beta equilibrium were also investigated by mode 
filling in Ref. [B] where the authors checked possible strong and weak hadronic decays (also 
multiple hadron decays) and found that slets stable against strong decays were most likely 
highly negatively charged. Mode filling studies are very important to show shell effect |38j, 
but it is difficult for large baryon numbers. Therefore, a multi-expansion liquid-drop model 
was developed [28J. Similar studies were done by He et al. at finite temperature [5 J and 
Zhang et al. with their suggested quark mass scaling (QMDTD) [39]. In this section, we 
apply our newly derived quark mass scaling in Eq. (j4]) to investigate the properties of slets. 

In relativistic heavy ion experiments, a slet, if formed, has no time to be in perfect beta 
equilibrium. We therefore regard it as a mixture of u, d and s quarks. Given three conserved 
quantities, i.e., the baryon number A, strangeness S and electric charge Z, we have the 
following equations: 

A = 1 -{N U + N d + N s ), (12) 
Z = 2 -N u - iiV d - ~iV s , (13) 
S = N s , (14) 

where N u , N d , N s are the number of u, d, and s quarks. For the mechanically stable slets, 
the internal pressure must be zero, i.e. 

P = 0. (15) 

We define the charge to baryon number ratio and strangeness fraction as f z = Z/A 
and f s = S/A = N s /A. Two different linear combinations of the Eqs. (fl2|) and (TIB"]) give 
N u = A + Z and N d + N s = 2A - Z. We then easily get 

U = NJA-1, (16) 
7 



fs 



3-{N u + N d )/A. 



Because we consider only < N u < 3A and < N d < 3A, the possible range for f z and 
f s are — 1 < f z < 2 and —3 < f s < 3. With a view to the relation f z + f s = 2 — N d /A, we 
have — 1 — f z < f s < 2 — f z if f z is fixed, and we can write — 1 — f s < f z < 2 — f s if f z is 
given. 
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Figure 1: The mechanically stable radii of 
a strangelet with baryon number A = 10 at 
temperature T = and 50 MeV is shown 
as functions of its strangeness fraction f s for 
given ratios of charge to baryon number. A 
minimum is reached on each curve. 
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Figure 2: The radius of a strangelet as func- 
tions of the charge to baryon number ratio at 
different given strangeness fraction. The up 
part is for A = 50 and T = 0, while the down 
part is for A = 5 and T = 50 MeV. One can 
also find a minimum on every curve. 



For a definite set of values for A, f z ,f s , one can easily get the quark numbers by N u = 
A(l + f z ), N d = A{2 — f z — f s ), N s = Af s . If we give an arbitrary value to the radius R, 
then the slet volume is V = (4/3)ttR 3 , and the density is n h = A/V = 3A/(47cR 3 ). At a 
given temperature, the quark density and temperature dependent mass can be accordingly 
calculated from Eq. (jl]). The respective chemical potential /x, {i = u,d,s) can then be 
obtained by solving rii = 3A^/(47ri? 3 ) with the expression in Eq. ([5]). The pressure and 
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free energy density are calculated by Eqs. (J6j) and (jHJ), respectively. We numerically vary R 
until the pressure becomes zero when the mechanically stable radius is reached. 

At the fixed baryon number A = 10, we give the slet radius as a function of the strangeness 
fraction f s at T = (the down part) and T = 50 MeV (the up part) in Fig. [U The solid, 
dashed, and dotted curves correspond, respectively, to the charge to baryon number ratio 
f z = —1, —0.5, and 0.2. It is obvious that the radius is not a monotonic function of 
strangeness. The position of the minimum radius depends on electric charge fraction. It is 
also natural that the radius increases with increasing temperature. 

Similarly in Fig. [2j we show the radius of strangelets as functions of the charge fraction. 
The variation of the radius with respect to electric charge is not monotonic at the fixed 
strangeness. The curves are also parabolas as in the above Fig. [TJ Generally, the minimum 
radius will appear at the middle position of the region (—1, 2 — /„). 

For comparison, we have also done calculations with the conventional bag model, with the 
results shown in Fig. [3j The left (right) panel shows the radius with respect to strangeness 
(charge) fraction for A = 10 at zero temperature. One also finds a minimum on each curve. 
Therefore, it is safe to conclude that the radius of slets is not monotonic. Instead, it has a 
minimum with respect to the charge and strangeness. This minimum is the smallest radius 
for a given pair of baryon number and temperature. However, the smallest radius does not 
necessarily correspond to the thermodynamically stable slets obtained by minimizing the 
free energy. 

Now we investigate the thermodynamically stable radius of a slet with given baryon num- 
ber at fixed temperature. In order to understand how the stable radius appears analytically, 
we may fall back on the fundamental differentiation equality of thermodynamics, i.e., 



dF = SdT - PdV + Y, VidNi. 




At fixed temperature T and P = 0, we have 



dF = Y,VidNi. 



(19) 
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Figure 3: The strangeness (left panel) and charge (right panel) dependence of the strangelet 
radius in the conventional bag model. The relevant parameters are A = 10, T = and 
5 1 / 4 = 180 MeV. A minimum is also obviously reached on every curve. 



If A and Z are fixed, we obtain naturally, 

d(F/A) 

Similarly, with fixed A and N s we get 

d(F/A) 



(20) 



(21) 



From Eqs. (12"U|) and (|2"T|) . we draw an important conclusion that the minimum of the free 
energy per baryon occur when = fi s with fixed A and Z in the F/A-f s panel, and when 
fJ'u = with fixed A and S in the F/A-f z panel. Therefore, we can get the stablest radius 
for a slet with given A and T by requiring the condition 



— — fJ'S 



(22) 



i.e., in this case only one chemical potential is independent. The only independent one can 
be determined by solving the equation 



(n u + n d + n s 



n h with n h 
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Figure 4: The free energy of strangelets 
versus the ratio of charge to baryon num- 
ber (left panel) and the strangeness fraction 
(right panel) for the two cases: A = 5, T = 
and A = 10, T = 50. 



Figure 5: Comparison of the smallest radius 
(the square) and the thermodynamically sta- 
ble radius (the triangle) in the R-f z plane 
(left panel) and in the R-f s (right panel) for 
A = 10 and T = 50 MeV. 



for an arbitrary radius. We finally vary the radius R so that the zero pressure condition is 
satisfied, and accordingly we obtain the stable radius. 

As mentioned in the introduction section, the size of slets is very important for analyzing 
their propagation and detection. Therefore, let's pay special attention to the case of zero 
temperature and try to derive an approximate expression for the radius. 



Because = J /if — m], Eq. (|22j) is equivalent to 



v l + m l = V v l + m\ = Ju 2 s + m 2 s = /I* 



(24) 



where we have used fi* as the common effective chemical potential. This means U{ = 
\J H* 2 — m 2 , or Xi = yj (/i*/mj) 2 — 1 where i = u, d, or s quarks. Substituting these into 
the Ui given in Eq. ((Sj) or Eq. at zero temperature, then substituting into Eq. (|23|) . we 
obtain an equation which contains the common chemical potential /i* and the radius R. Sim- 
ilarly substitution into the pressure expression in Eq. §6§ or (TTTj) . we get another equation 
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of //* and R. The radius is then obtained by solving the two equations of /i* and R. 

Because the surface and curvature terms can be regarded as a perturbation to the volume 
terms, and also because the interaction part of the quark mass is greater than the quark 
current masses (mi > m^o), we can ignore the finite-size effect and iso-spin effect to derive 
a first-order approximation for the radius. In this case, we have v u = v& = v s = u and 
x u = Xd = x s = xq, and the zero pressure condition becomes 



x (2xq — 3)Jxq + 1 + 3arcsh(x ) — I2z x Jxq + 1 — arcsh(x ) 



(25) 



The positive solution of this equation is Xq ~ 1.3278478 (The trivial solution x = and 
the non-physical solution xq = —1.3278478 were discarded). At the same time, Eqs. ([9]) and 
( 1231) at the same order approximation gives Hi = z^f/vr 2 = = 3A/(4irR 3 ). Combining this 
with x = v^/mi = (7r 2 nb)^ 3 / (D /n]/ 3 ) = (7m b ) 2 / 3 /D, or n h = (x -D) 3//2 /vr, we immediately 
have 

R = W^-A^ = r slet A^. (26) 
For the chosen value D = (156 MeV) 2 , the reduced slet radius is 

rsiet = R/A 1/3 = (3/4) 1 / 3 / yx^D = 0.9973 fm. (27) 

This value is smaller than that of normal nuclei, but bigger than the recent value (~ 0.94 
fm) [40J from the conventional bag model calculations. Equation ( 1261) shows that the slet 
radius is inversely proportional to the square root of the confinement parameter D. If one 
uses a bigger D value, then the slet radius can be very small. In that case, however, SQM 
will not be absolutely stable. 

It should be emphasized that the expression for the slet radius in Eq. (1261) is only the 
lowest order approximation because it was obtained by ignoring the finite size effect and 
isospin dependence (quark mass difference). The actual size is a a little bit bigger, from 1 
fm (for large baryon numbers) to about 1.2 fm (for small baryon numbers). So Eq. (1261) is 
only accurate for slets with large baryon numbers. In the following we continue to present 
the numerical results. 
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Figure 6: The baryon number dependence of the strangelet radius at temperature T = 
(solid curve), 20 MeV (dashed curve), and 50 MeV (dotted curve). The dash-dotted curve 
is for the ordinary nuclei radius R = 1.12A 1//3 fm. 

We show all possible slets with positive strangeness in Fig. HI The full dots are for A = 5 
and T = 0, while the open circles are for A = 10 and T = 50 MeV. The left (right) panel 
gives the distribution of the free energy per baryon on the f z (f s ) axis. The triangle denotes 
the minimum of the free energy per baryon, i.e., the most stable slets. For a given pair 
of A and T, therefore, we have two special values for the radius of a slet via varying the 
strangeness fraction f s and the charge to baryon number ratio f z : the first one is obtained 
by minimizing the radius itself, and thus is the smallest radius, while the second one is 
obtained by minimizing the free energy of the system, and is thus the thermodynamically 
stable radius. The obvious difference between the smallest radius (square) and the stable 
radius (triangle) is compared in Fig. [5l where the smallest radius is 2.95 fm for A = 10 and 
T = 50 MeV while the stable radius is 3.05 fm. 

The baryon number dependence of the stable slet radius is shown in Fig. [6] at temperature 
T = (solid line), 20 MeV (dashed line), and 50 MeV (dotted line). To compare, the normal 
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nuclear radii have also been shown with a dash-dotted line. We see that the radii of slets 
are comparable with those of ordinary nuclei. 
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Figure 7: The free energy per baryon (left 
axis) and the charge to baryon number Z/A 
(right axis) of stable strangelets vary with 
the baryon number A. They are generally 
decreasing functions. With increasing tem- 
perature, the decreasing speed goes down, 
and Z/A is very flat at high temperature. 
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Figure 8: The strangeness per baryon for 
stable strangelets as a function of the baryon 
number at zero temperature. For very small 
baryon numbers, the strangeness fraction in 
the present model (solid curve) becomes neg- 
ative, while that in the bag model (dotted 
curve) is always positive. 



Fig. [7] gives the free energy per baryon and the charge to baryon number ratio of stable 
slets as a function of the baryon number. It is seen that the free energy per baryon decreases 
with both increasing temperature and baryon number. When the temperature is as high as 
up to 50 MeV, the dash-dotted line for charge to baryon number becomes very flat. The 
smaller the baryon number, the more important the influence of temperature on the charge 
to baryon number. 

In Fig. [HI we show the strangeness fraction f s , with a solid line, as a function of the 
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baryon number at zero temperature. The striking feature is that the strangeness fraction 
becomes negative for very small baryon numbers. In fact, this feature has already been 
seen in the up part of the right panel in Fig. H] where the minimum of the free energy was 
not reached at positive strangeness. For comparison, we also give, in the same figure with 
a dotted line, the results from the conventional bag model in which the quark masses are 
constant and the strangeness fraction is always positive. In the present model, however, the 
quark masses are density and temperature dependent. Therefore, the strangeness fraction 
becomes negative for very small baryon numbers, and the charge to baryon number ratio 
becomes bigger than that in the bag model. 



4 Summary 

We have studied the properties of strangelets in a new quark mass scaling without imposing 
beta equilibrium. It is found that the radius of strangelets is not a monotonic increasing or 
decreasing function of either electric charge or strangeness. By varying the strangeness and 
charge for a given pair of baryon number and temperature, we have calculated the smallest ra- 
dius and thermodynamically stable radius, and have shown that they are generally different. 
The stable radii of strangelets can be calculated approximately by the law R = r s \ ct A 1 ^. If 
SQM is absolutely stable, the reduced radius of strangelets is r slet = (3/4) 1 / 3 / \/xqD ~ 1 fm. 
The smallest radius appears always at positive strangeness. However, the stable strangelets 
could have negative strangeness fraction for very small baryon numbers. This means that 
very small strangelets may contain anti-strangeness. 

The present study is still model-dependent. Many problems, e.g., how the anti-strangeness 
content influences the stability and detectability of strangelets in heavy ion experiments, to 
what extent the behavior of the radius can be used to analyze the production and detection 
of strangelets, etc, need to be further explored. 
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